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Abstract. In order to specify cosmologically motivated initial conditions for major galaxy mergers (mass ratios < 4: 1) that are 
supposed to explain the formation of elliptical galaxies we study the orbital parameters of major mergers of cold dark matter 
halos using a high-resolution cosmological simulation. Almost half of all encounters are nearly parabolic with eccentricities 
e ~ 1 and no correlations between the halo spin planes or the orbital planes. The pericentric argument CO shows no correlation 
with the other orbital parameters and is distributed randomly. In addition we find that 50 % of typical pericenter distances are 
larger than half the halo's virial radii which is much larger than typically assumed in numerical simulations of galaxy mergers. 
In contrast to the usual assumption made in semi-analytic models of galaxy formation the circularities of major mergers are 
found to be not randomly distributed but to peak around a value of e ~ 0.5. Additionally all results are independent of the 
minimum progenitor mass and major merger definitions (i.e. mass ratios < 4:1; 3:1; 2:1). 
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1. Introduction 



iToomre & Toomrj ( 1 19721) (TT72) proposed that elliptical 
galaxies could form from mergers of disk galaxies. Thereafter 
this "merger hypothesis" has been investigated in great de- 
tails by many authors (for recent reviews see Burkert & Naab 
2003a,b) using numerical simulations. The first fully self- 
consistent models of two equal-mass stellar disks embedded 
in dark matter haloswere p erformed by jN egroponte & White 
(Il983h . iBarned (Il988l) and iHernauistI ( Il992l) . Typically, the 
spiral galaxies were constructed in dynamical equilibrium 
(Hernauist, 1993a) and approached each other on parabolic or- 
bits with pericenter distances of twice the scale length of the 
larger disk galaxy. Already the first simulations demonstrated 
that mergers lead to slowly rotating, pressure supported and 
anisotropic systems that generally follow a de Vauco uleurs sur- 
face density profile and resemble elliptical galaxies. iHernauisl 
( 1993b) noted that a massive central bulge component was re- 
quired in the progenitors in order to explain the high stellar 
densities observed in the cores of ellipticals. More recent in- 
vestigations showed that the initial mass ratio of the merging 
spirals determines the isophotal shapes and kine matical prop- 
erties of the merger remnants ( Bendo & Earned, EoOO ) and it 
was suggested e.g. by Naab. Burkert. & Hernauist ( 1999) that 
the initial mass ratio might explain the observed dichotomy be- 
tween disky, fast rotating and boxy, slowly rotating elliptic als 
(Bender, Doebereiner, & Moellenhoff, 1988; BendeJ. ll98^ . 

Recently, iNaab & BurkertI ( l2003l) performed the to date 
largest parameter survey of collisionless binary disk mergers 
with varying mass ratios, different initial disk orientations and 



orbital parameters. They found that the fine structure of the 
merger remnant in general depends strongly on the adopted ini- 
tial conditions. Certain initial conditions in fact lead to disky, 
anisotropic remnants which are not observed. Clearly, a deeper 
understanding of the origin of the early type family of galaxies 
not only requires a more sophisticated treatment of the various 
physical processes but also a better understanding of the initial 
orbital configurations of the interacting galaxies as predicted 
by currently favoured cosmological models. 

Up to now due to the lack of knowledge of the appropri- 
ate initial conditions, simulations were forced to cover a wide 
parameter space by setting up mergers with different orbital 
configurations. This approach also has the drawback that it is 
not clear how relevant a given parameter combination is which 
leads to a certain fine structure of the merger remnants. In 
this paper we therefore present a detailed analysis of a high- 
resolution, large-scale cosmological simulation that was car- 
ried out by the Virgo Supercomputing Consortium using com- 
puters based at Computing Centre of the Max-Planck Society 
in Garching and at the Edinburgh Parallel Computing Centre,' 
and derive the typical orbital parameters predicted by cosmol- 
ogy for merging dark matter halos which should also be repre- 
sentative for their visible galactic components. 

1.1. The cosmological simulation 

The simulation we analy se wa s performed within the GIF 
project (Kauffmann et al l ll999l) 



and we will here focus on 
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the ACDM simulation. This simulation was carried out in a 
box of size 141.3 Mpc with 256^* particles each hav- 
ing a mass of 1.4 x IO'^Mq/z^' and cosmological parameters 
Q.A — 0.7, Q.() — 0.3, ag = 0.9, h — 0.7 and a force softening of 
20 Kpc//i. The softening leads to an error of a few percent in 
the calculation of the orbital parameters. The positions and ve- 
locities of the particles have been saved at 44 different red- 
shifts. Additionally, for each output redshift, a list of halo prop- 
erties have been provided by the GIF project. When two ha- 
los approach each other their orbit is going to change due to 
the transfer of orbital angular momentum to the halo's internal 
angular momentum which in the following is called spin and 
should not be confused with the cosmological spin parameter 
defined e.g. in Peebles ( 1993). The determination of the initial 
orbital parameters therefore is a question of the 'right timing'. 
We try to identify these parameters at a time when the inter- 
action between the halos is still weak and one can assume a 
Keplerian two-body situation, using the positions of the most 
bound particles of each halo. At each redshift we go through the 
list of halos identified by using the friends-of-friends (FOF) al- 
gorithm and identify the positions of the most bound particles. 
If at one redshift a halo has disappeared through merging with 
another halo, we look up the position of its most bound par- 
ticle at the previous redshift and check whether the distance 
to the most bound particle of the other halo, with which it is 
going to merge, is larger than the sum of both virial radii. If 
so, we derive the orbital informations using the data from this 
redshift. To make sure that the merger is not just a flyby of 
an unbound halo we check at a redshift later than the merging 
redshift whether the separation of the most bound particles has 
increased again. Since outputs are generated roughly every 0.5 
Gyr we expect not to miss even very eccentric orbits. 

2. The reduced two-body problem 

The orbital parameters are defined as specified in the pioneer- 
ing work of TT72. We simplify the problem by reducing the 
two-halo system to a two-body system with each halo being 
represented by its most bound particle which sits in the poten- 
tial minimum and is expected to be the most 'stable' particle, 
allowing us to follow the orbital evolution of the progenitor 
halos during the early stages of the merger even when a clear 
definition of the centre of mass of the extended and already 
partly overlapping progenitor halos is not possible anymore. 
We checked that the position of the most bound particle is not 
very different from the centre of mass of each halo when the 
halos are well separated and that it agrees with the region of 
highest dark matter density. We additionally tested if the re- 
sults presented in this paper show any significant dependence 
on using the centre of mass of each halo, calculated by averag- 
ing the particle positions inside different fractions of the virial 
radius, instead of the position of the most bound particle and 
found no significant difference. The position of the most bound 
particle in each of the halos will be denoted by r with the index 
h, like halo, for the more massive of the two merging halos and 
index s, like satellite, for the less massive halo. Throughout this 
paper we will use these indexes to refer quantities to the more 
and less massive of the merging halos. 



In the following we will assume that each of the two point 
masses moves with the average velocity of the particles inside 
the virial radius of the corresponding halo. 

Knowing the average velocity and mass inside the virial 
radii of the merging halos it is possible to solve the gravita- 
tional two body problem analytically by making the transition 
to the equivalent one -body probl em of the reduced mass system 
(e.g. Goldstein, Pool e, & Safk3l2002l) . The equation of motion 
reads 



AT 



,M/,M.v r 



(1) 



with M/, and being the mass inside the virial radius of the 
more and less massive of the merging halos, respectively. The 
reduced mass is defined by 

MhMs 



(2) 



M,,+M, 

Equation^describes the motion of a fictitious particle of mass 
jj moving around a point at a distance r = |r/, — r^l- To spec- 
ify the orbit we calculate the total energy and orbital angular 
momentum according to 



P 1 
E ^ -^r 

and 



L = jjrxr, 



2 GM,M 



(3) 



(4) 



respectively. In this study we also investigate the pericenter dis- 
tance rperi and the eccentricity e of orbits, which we can readily 
calculate from the energy and orbital angular momentum by 



2£L2 



and 



' pen 



(l+e)ijGMiM 



(5) 



(6) 



In the following we will identify orbits according to their to- 
tal energy and eccentricity into hyperbolic {E > 0, e > 1), 
paraboUc (£ = 0, e = 1) and elliptic (E <0, e < 1) orbits. 

2.1. Limits on tine two-body approximation 

As summarised in the previous section we assume the two 
merging halos to be fairly isolated and without experiencing 
any other severe interactions before they merge. This however, 
is not guaranteed and needs to be checked. We calculated the 
total energy and the orbital angular momentum of the orbits of 
the two halos at two simulation outputs before they merge. To 
do so we followed the most massive progenitors backward in 
time and analysed their orbital parameters, including all halos 
with more than 50 particles and with a mass ratio of 4: 1 or less. 
The filled circles in top left panel of Fig.^show the total energy 
of the merging halos E„ and E„^i, one and two time steps be- 
fore the merger, respectively. As can be seen the energies seem 
to be not too different between the two time steps, with a ten- 
dency for bound orbits at time step « — 1 to become more bound 
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at n. The dashed histogram in top right panel shows the distri- 
bution of ratios E„^i/E„ for all identified mergers. The dis- 
tribution is skewed toward values smaller than one which con- 
firms the tendency of the orbits to become more bound between 
two time steps. The filled circles in the left middle panel of the 
same figure shows the correlation between the orbital angular 
momentum at two preceding time steps before the merger The 
dashed histogram in the right middle panel shows the distribu- 
tion of the ratio of these orbital angular momenta. In general 
the orbital angular momentum between two outputs does not 
change. However we find that there is a considerable scatter 
which is symmetrically distributed around L„ = L„_i. 

The aim of this work is to present initial conditions for 
high resolution simulations of isolated merging galaxies. This 
means that we ideally would like the energy and angular mo- 
mentum to be conserved or not to change significantly while 
the two merging halos are well separated on their approach. As 
can be seen in the lower left panel of Fig.^many of the pro- 
genitors are not isolated in the sense that they do not accrete 
further mass before they merge. To choose fairly isolated halos 
we therefore decide to only consider those mergers in which 
l.ll(M/,,„ +M,,„) > Mh,n-i +M,,„_i > 0.9(M/,,„ +M,,„). The 
result of this selection criterion is shown by the grey triangles 
in the left panels of Fig. [J In all cases the correlations between 
the quantities at the different redshifts become tighter. In ad- 
dition we find that the distribution in the ratio of energies is 
now centred around a value of one (solid histogram right top 
panel). The shape of the distribution in the ratios of the angu- 
lar momenta does not change much though it becomes a bit 
narrower. In the following we will use above mass criterion to 
select mergers for the sample used here. This criteria roughly 
rejects 50% of the mergers detected in the simulation. 

3. Orbital parameters rpen and e for major mergers 

We start by analysing the dependencies and correlations of 
the pericenter distance rperi and the eccentricity e of merging 
halos extracted from the simulation data at one output prior 
to the merger We restrict ourselves to mergers with mass ra- 
tios of 4: 1 or less as numerical models have shown that larger 
mass ratios do not resemble ellipti cal galaxies anymore as the 
more massive disk is not destroyed jBarnes & Hernquislll992t 
lNaab& Burkert .200311) . 

In figure |2] the distribution of eccentricities e of merging 
halos as function of the minimum mass of the progenitor ha- 
los is shown. Roughly half of the orbits are parabolic or very 
close to parabolic. We find ~ 40% of the orbits in the range 
e = 1 ±0.1. This result is independent of the minimum mass 
cut applied. N-body simulations of merging galaxies assume 
parabolic orbits (e.g. iBarnesL Il988l) . which is only the case 
for half of the mergers we find. In fig. |3lwe show the depen- 
dence of our results on the mass ratio of the major mergers. 
We find no dependence with the mass ratio, which indicates a 
self-similarity of the formation process. As many characteristic 
properties of elliptical galaxies as e.g. the amount of rotation 
are determined by the mass ratio (Naab & Burkert 2003) the 
cosmological models would predict that massive and low-mass 
ellipticals should have similar properties which is not in agree- 




pail,n-l parL,n-l part,n 



Fig. 1. Upper left panel: Total energy of the two halo sys- 
tem. Upper right panel: Distribution of ratios E„- 1 /E„. Middle 
panel: Orbital angular momentum of the two halo system. 
Middle right panel: Distribution of ratios L„_i/L„. Lower 
panel: Mass of the progenitors in units of particle mass 
in the simulation. Lower right panel: Distribution of ratios 
N part, n-\/N part. n- The indcxcs n and n — 1 refer to one and two 
simulation time steps before the merger. The circles show all 
identified mergers, and the triangles only those when applying 
the mass criterion. Dashed histograms show the distribution for 
all mergers and solid histograms only for those satisfying the 
mass criterion. 



ment with observations. In contrast with the theoretical pre- 
dictions luminous ellipticals are preferentially anisotropic and 
slowly rotating while less luminous ellipticals rotate fast and 
appear isotropic (Bender 1988). The distribution of the peri- 
center distance in units of the virial radius of the more massive 
progenitor halo and its dependence on minimum mass is shown 
in the left panels of fig.|3 The distribution shows no significant 
dependence on the minimum mass. Small pericenter are more 
frequent than larger ones. This is actually what one would ex- 
pect, because halos which are on orbits leading to a very close 
encounter are more likely to merge than those which pass each 
other from very far. Roughly 85 % of the mergers have peri- 
center distances which are larger than 0.1 ry,,^/,. Fig. EJshow 




Khochfar & Burkert: Merging Dark Matter Halos 

0.25 



.2 

?3 



0.2 
j0.15 
I 0.1 
0.05 

0.2 
1.15 
0.1 
0.05 

0.2 
;0.15 
I 0.1 
0.05 




M, /M < 4 

h s 




M, /M < 3 

h s 





0.5 1 1.5 2 




0.5 1 1.5 



Fig. 2. Left panel: Eccentricity distribution of merging halos Fig. 3. Left panel: Fraction of merging halos with mass larger 
with mass larger than Npart ■ 2 x lO'^Mpj. Npan is the number of than lO'^M© and different mass rations M^/Ms, on initial or- 
dark matter particles. Right panel: Corresponding cumulative bits with eccentricity e. Right panel: Corresponding cumulative 



fraction of eccentricities. 



fraction of eccentricities. 



that the dependence on the mass ratio of the mergers is neg- 
ligibly small. Merger simulations usually adopt initial condi- 
tions with small pericenter distances which are typically of or- 
der 0.05r,,,>./i, leading to fast merging that saves computational 
time. This assumption is not in agreement with cosmological 
models and might affect the inferred structure of merger rem- 
nants using numerical simulations. According to the reduced 
two-body problem rperi °^ L^, which means that the orbital an- 
gular momentum in typical merger simulations is less than that 
for merging halos in self-consistent cosmological simulations. 
Less angular momentum will therefore be transferred during 
the merger process and the structure of the remnant could be 
different. In fig.|5]we correlate the pericenter distance in units 
of Mpc to the pericenter distance in units of r,,,>/,. We find that 
on average roughly 50% of all mergers have pericenter dis- 
tances less than 0.5r,„./,. About 70% of all mergers lead to a 
first pericenter passage which penetrates the virial radius of the 
larger halo r,,,,./,. 

Since we have found no dependence on the minimum mass 
and on the mass ratio we continue our investigations using as 
a standard assumption M„„„ = lO'^M© which corresponds to 
the typical halo size of massive galaxies and M^/Ms < 4. Not 
every random orbit is going to lead to a merger and it is im- 



portant to see if a preferred orbit configuration exists leading 
to mergers. In fig. |6l we correlate rp^i-i and e. Mergers with 
rperi < O-ir^irji are almost all on parabolic orbits with e ^ I. 
Orbits with rperi > 0.1r,,,>7! have a scatter in e which increases 
with pericenter distance toward bound orbits. The same be- 
haviour is found looking at the correlation between eccentrici- 
ties and pericenter distances in units of Mpc. 

The expectation from the two-body solution is that the spe- 
cific angular momentum L^p = L//j for remnants of the same 

1/2 

mass Mil + Ms is proportional to r^J^^- with a scatter because of 
different eccentricities of the orbits. In fig.0we show the cor- 
relations found between these two quantities. The lines in the 
left and right panel of fig.^lare power law fits to the data with 



L,p = 1.11 X 10 



-10 



' pen 
^virji 



0.48 



Mpc^ 
yr 



for the pericenter distance in units of r,,,>,/j and 

0-53 i^_„2 



L,p = 2.84 X 10 



10 / 'pen 



' Mpc^ 



(7) 



(8) 



Mpc / yr 

for the pericenter in units of Mpc. The fits show that the data 

1/2 

is following the trend of Lyp °^ r '■. The larger scatter in the 



Khochfar & Burkert: Merging Dark Matter Halos 



5 




r . / r . , 

pen vir,h 



Fig. 4. Left panel: Fraction of merging halos with different 
pericenter distances in units of the more massive progeni- 
tor's virial radius r^irji and with progenitor masses larger than 
Npart '2 X IO^'Mq. The corresponding cumulative fraction of 
pericenter distances are also shown. All results are shown for 
mergers with a mass ratio of M/,/Ms < 4. Right panel: Fraction 
of merging halos with mass larger than lO^^M© and differ- 
ent mass ratios Mi, /Ms, on orbits with different pericenter dis- 
tances in units of the more massive progenitor's virial radius 



correlation with rperi in units of r,,,>/j is due to additional scatter 
introduced on the horizontal axis by the variation in the values 
of ryir,h at a given 

^ pen- 
Hyperbolic encounters are often characterised by the im- 
pact parameter b. The impact parameter is defined by the spa- 
tial separation perpendicular to the initial velocity V(t = —°°) 
which has the value 



(9) 



Because of the conservation of orbital angular momentum the 
impact parameter is given by 



(10) 
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Fig. 5. Pericenter distance in units of Mpc against pericen- 
ter distance in units of /•,,,>,/;■ The values at each line indicate 
the completeness limit, meaning the fraction of mergers with 
fperi/fvir.h < 0.5 and 1. Results are shown for progenitor ha- 
los of mass larger than IO^^Mq and mergers with mass ratio 
Ml, /Ms < 4. 
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Fig. 6. Correlation between the eccentricity and pericenter dis- 
tance of merging halos found in the simulation. Results are 
shown for progenitor halos of mass larger than lO^^M© and 
mergers with mass ratio Mi, /Ms < 4. 



Figure |S] shows the probability distribution function of impact 
parameters which shows a peak around b w 0.5. We fit the dis- 
tribution with a function of the form: 

"0-1 



exp 



dx, 



(11) 



dP J 1 flo 

dx fl2 fli \fli y ' \ ai 

and we present the best fit parameters in tabled In the upper 
graph of fig.|9]we show the correlation between the impact pa- 
rameter and the pericenter distance. The line is a power-law fit 
to the data with 



' pen 



- 0.3b' 



0.72 



(12) 



The data shows a large scatter which makes the power-law fit a 
rather poor fit. However, to guide the eye and show the general 
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Fig. 7. Correlation between the specific angular momentum Lsp 
of the orbit and its pericenter distance in units of rv,>,/i (left 
panel) and Mpc (right panel). Lines represent power law fits to 
the data. Results are shown for progenitor halos of mass larger 
than lO'^Mg, and mergers with mass ratio M/,/Ms < 4. 
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Fig. 8. Upper graph: Probability density of hyperbolic orbits 
with different impact parameter b. Lower graph: To the upper 
graph corresponding cumulative fraction. Results are shown for 
progenitor halos of mass larger than lO'^Mg, and mergers with 
mass ratio MjjMs < 4. 



trend we include this fit. The lower graph of fig. |9] displays 
the correlation between eccentricity and impact parameter of 
the encounter Again it becomes evident that the many of the 
orbits are close to be parabolic and that only a small fraction 
is significant different from parabolic. Encounters having 
e > 1 merge very slowly if at all. This explains why only 
hyperbolic encounters with small impact parameters b < I 
Mpc, corresponding to a close flyby, lead to mergers. 

Another parameter commonly used to describe bound or- 
bits with £■ < is the circularity e defined as the ratio of the or- 
bital angular momentum to the angular momentum of a circular 
orbit with the same energy. The circularity of a bound orbit can 
be derived applying the virial theorem U = —2T. With 



GMiM 



Vn 



-2E 



' arc — T J ' arc 

2E 



and the angular momentum 

^circ — ^ circt^circ : 

the circularity becomes 
^_ L L 



(13) 



(14) 



(15) 



From eq. E]one sees that the circularity can only be defined 
for orbits with E <0. Manipulating the appropriate equations 
gives the following relation for the circularity and eccentricity 
of an orbit: 



(16) 



The upper graph of fig. [TOl presents the distri bution of circu- 
larities found. This result agrees with that of iTormenl ( 1 19971) 
who found the circularities to be distributed with a peak around 
e « 0.5. However, in their analyses they considered only minor 
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Fig. 9. Upper graph: Correlation between pericenter distance 
and impact parameter of hyperbolic orbits leading to merg- 
ers. Lower graph: Correlation between the eccentricity and im- 
pact parameter corresponding to the orbits in the upper graph. 
Results are shown for progenitor halos of mass larger than 
IO^^Mq and mergers with mass ratio MjjMs < 4. 



mergers Mj,/Ms ^ 4 in a cluster environment, where the en- 
hanced gravitational field might lead to a change of the circular- 
ity distribution. Our results show that these effects seem to play 
not an important role, and that the circularity distribution seems 
to be universal. An important consequence for semi-analytic 
modelling is that the circularity e used to calculate the dynam- 
ical friction time scale for major mergers m ust not be drawn 
from a uniform distribution as e.g. done in iKauffmann et all 
(il999i) but from the distribution found here. 
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Fig. 10. Upper graph: Probability density of bound orbits with 
different circularity. Lower graph: To the upper graph corre- 
sponding cumulative fraction. Results are shown for progenitor 
halos of mass larger than lO'^M© and mergers with mass ratio 
Mh/M, < 4. 
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Fig. 11. Correlation between pericenter distance and circular- 
ity of bound orbits leading to mergers. Results are shown for 
progenitor halos of mass larger than 4 x lO'^M© and mergers 
with mass ration Mi, /Ms < 4. 



Fig- E] shows that the correlation between circularity and 
the pericenter distance. We find that the data can be fitted well 
by following power-law 

=0.38e2'^Mpc (17) 
4. Parameters co & / 

If halos have spin, the orbital parameters discussed in the last 
section are not enough to fully describe the geometry of the 
encounter. Additional constraints on the position of the spin 
vectors S of both halos are required. These spin vectors are in 
general expected to align with the spin vectors of the galactic 
disks that form by gas infall and kinetic energy dissipation. In 
calculating the spin vectors we follow iLemson & Kauffmanni 



Orbit 



Spin ' 
Plane 




Fig. 12. Definition of the angles / and co following TT72 (figure 
reproduced with permission of A. Toomre). 



who derived spin parameters from a similar GIF simula- 
tion. They found that for halos with more than 70 particles halo 
properties are calculated reliably. Additionally they require that 
between 35 and 95 per cent of the friends-of-friends mass of the 
halo are inside of the virial radius. They find that 80 per cent of 
the halos satisfy these requirements. In our analysis we adopt 
limits of 50 particles and more per halo and between 50 and 95 
per cent of the friends-of-friends mass inside the virial radius. 
In that way we make sure that we get results which are not in- 
fluenced by ambiguities in the way we treat the particles inside 
each halo. In fig.^lthe definition of the two corresponding an- 
gles is shown. The angle / is defined in the rest frame of the 
halo as the angle between the spin plane of the halo and the or- 
bital plane of the satellite; and in the rest frame of the satellite 
as the angle between the spin plane of the satellite and the or- 
bital plane of the halo. By spin and orbital plane we refer to the 
planes perpendicular to the spin and orbital angular momentum 
vector, respectively. These two angles if, and are independent 
and by definition |/| < 180°, where ; = 0° is a prograde and 
i — 180° a retrograde encounter. Additionally, the pericentric 
argument CO, is defined as the angle between the line of nodes 
and separation vector at pericenter, and has values ranging from 
CO = -90° to CO = 90°. It is not defined for / = 0° or / = 180°. 

Figures [O]- [21 show the correlation between the angles ; 
and CO and the minimum mass of progenitor halos or the mass 
ratios of the mergers. The angles //, and is are distributed fol- 
lowing a sinus function. The fit gets naturally poorer at high 
minimum masses because of the smaller sample of halos. The 
solid lines in figure [T3landll4lare fits of the form oc \sin{x)\. If 
the angle between two vectors is sinus-distributed, the two vec- 
tors have no correlation. This can be understood from looking 
at the probability of drawing a random vector pointing from the 
centre of a sphere to its surface. If every point on the surface is 
equally likely, the probability of finding an angle / for example 
between the x-axis and a random vector will be proportional 
to sin(/). We therefore conclude that the spin plane and the or- 
bital angular momentum plane have no correlation. Inspection 
of figure [21i"eveals that the angle co is equally distributed and 
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Fig. 13. Left column: Angle between halo spin plan and orbital 
plane for different choices of minimum progenitor mass and a 
major merger definition of Mh/Ms < 4. Right column: Same as 
left column but now for the angle between the spin plane of the 
satellite and orbital plane. 



Fig. 14. Left column: Angle between halo spin plan and orbital 
plane for different choices of major merger definition and fixed 
minimum progenitor mass of lO'^M©. Right column: Same as 
left column but now for the angle between the spin plane of the 
satellite and orbital plane. 



independent of the minimum halo mass and major merger def- 
inition. In contrast to /, (O is defined by the position of a vector 
lying in a predefined plane. As a result an equal distribution in 
CO means that this angle is randomly distributed. 

Defining K as the angle between the spin planes, one would 
expect from the results presented above that the spin vectors 
are also not correlated with each other, and that K should be 
sinus-like distributed. The distribution of K is shown in fig. 1161 
and is indeed sinus like. 

It is important to investigate possible correlations between 
the orbital parameters r^g,, and e and the angles introduced 
above. Fig. [n] illustrates the correlation between the angle 
is and the orbital parameters. The orbital parameters are not 
correlated with /j. The same results are found for //,. The 
pericentric argument cOj shows no correlation with the orbital 
parameters as shown in fig. ^] The results are the same for 

(OA. 



Semi-analytic mo dels llMaller. Dekel. & SomervillS 120021: 
IVitvitska et al.[ l2002l) assume that during mergers the orbital 
angular momentum gets transformed into spin of the remnant 
halo. These models and models in which angular momentum 
is acquired by tidal torques e.g. Porciani, Dekel & Hoffmann 
(2002) reproduce the spin distribution of halos found in N-body 
simulations. However the merger picture for the build up of 
halo spins uses some assumptions which still need confirmation 
by N-body simulations. As a first step the amount of angular 
momentum in the orbit must be investigated. Fig.^]shows the 
distribution of the fraction of orbital angular momentum to spin 
of the halos Sh and spin of the satellite Ss. 

iMaller. Dek el. & Somerville (200Z) define a parameter / 
for mergers 



VyirMrvir.h/J 

with Viirji as the circular velocity of the more massive progen- 
itor. The value of this parameter is set to be / ~ 0.42 for their 
model, in which spin is acquired from orbital angular momen- 
tum. In Fig. |20lthe distribution of / is displayed. We fit the 
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Fig. 15. Left column: Angle between pericenter vector and 
node line for different choices of minimum progenitor mass 
and a major merger definition of M/j/M^ < 4. Right column: 
Same as left column but now for fixed minimum progenitor 



mass of lO^^Mp, and different choices of major merger defini- 
tion. Results are the same for CO; and CO/,. 




Fig. 17. Upper graph: Correlation between of the satel- 
lite and the eccentricities of orbits leading to mergers. Lower 
graph: Correlation between the rperi and is corresponding to 
the orbits in the upper graph. Results are shown for progenitor 
halos of mass larger than 10 '^M© and mergers with mass ratio 
Mh/Ms < 4. 
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Fig. 16. The distribution of angles K between the two spin 
planes. Results are shown for progenitor halos of mass larger 
than lO'^M© and mergers with mass ratio Mi,/Ms < 4. 



Fig. 18. Upper graph: Correlation between co.s and e for orbits 
leading to mergers. Lower graph: Correlation between rperi and 
cOi, corresponding to the orbits in the upper graph. Results are 
shown for progenitor halos of mass larger than IO^^Mq and 
mergers with mass ratio Mh/Ms < 4. 



distributions in Fig. ^]&|20| with functions of the form pre- 
sented in eq. ^\ and the fitting parameters can be found in 
tabled The distribution for / peaks at 0.6 and has its mean 
at ~ 1 . which disagrees with the value required in the model of 
iMaller. Dekel. & Somervill&c200a) . 

5. Discussion and Conclusions 

We have investigated the orbital parameters of merging dark 
matter halos in a cosmological large scale simulation finding 
that the parameters used generally in simulations of merging 
galaxies are in agreement with those found in the large scale 
simulation. However, parameters used in merger simulations up 
to now only occupy a very small part of the region of parame- 
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X 


ao 


ai 


02 


b 


2.636 


0.753 


0.326 


L/Sh 


1.727 


7.408 


0.662 


L/S, 


1.624 


14.132 


0.726 


f 


2.638 


1.2 


0.307 



Table 1. Best-fit parameters using eq. ^Jto the distributions 
shown inFig.|8l[ll&|20| 
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Fig. 19. The upper two graphs show the distribution and cor- 
responding cumulative fraction of mergers with different L/Si,. 
The solid line in the upper of the two graphs is the fit using 
eq. ^2 Lower two graphs show the same as the upper graphs 
but now for the fraction L/Ss- Results are shown for progenitor 
halos of mass larger than IO^^Mq and mergers with mass ratio 
M,,/M, < 4. 



ter space covered by the results from the large scale simulation. 
In particular in most cases the calculated pericenter distances 
are more than a factor two larger than those used generally in 
merger simulations. First results from simulations using larger 
pericenter distances indicate more isotropic remnants, which 
agree better with observations (Naab & Burkert, in prepara- 
tion). About half of the orbits are close to parabolic and those 
which are not show larger pericenter distances on average. The 
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Fig. 20. Upper graph: Distribution of parameter / found in the 
simulations. Lower graph: Corresponding cumulative fraction 
found in the simulation. Results are shown for progenitor ha- 



tight correlation rp 



2.17 



which we find between the peri- 



los of mass larger than 10 M r: and mergers with mass ration 
Mi,/M, < 4. 



center distance and the circularity just supports that we have 
chosen halos when they are still well separated and not much 
interaction has taken place. Interaction leads to a circulariza- 
tion of the orbit predicting an increase of the circularity e on 
average with decreasing pericenter distance rp^,., . On the other 
hand hyperbolic orbits need to have very high kinetic energies 
to allow for small pericenter distances which the simulation 
shows to be unlikely. The pericenter distances however, show 
in the case of hyperbolic orbits a correlation with the impact 
parameter of the form rperi °^ b^ ''^, allowing to disentangle the 
different contributions to the orbital angular momentum using 
the additional correlation L^p °^ r'^'^^j . 

The distribution of circularities found is in contrast to what 
is commonly assumed in semi-analytic models of galaxy for- 
mation, which assume a random distribution of circularities for 
the orbits of infalling satellites. Instead our results indicate a 
distribution with a peak around £ « 0.5. 

We investigated the correlation between the spin plane of 
the main halo and the spin plane of the satellite with the or- 
bital plane of the merger, finding no correlation. The planes 
are aligned randomly which gives statistical weight to certain 
merger geometries. Additionally the spin planes are not corre- 
lated with each other and the angle ; between the spin planes 
and the orbital plane shows no correlation with the other orbital 
parameters rperi and e which justifies its use as an independent 
initial parameter. The pericentric argument CO follows a random 
distribution and shows no correlation with the other orbital pa- 
rameters rperi and e. 

We te sted the requirement needed in th e semi-analytic 
model of (Mailer, Dekel, & Somerville, '2002*) for the acquir- 
ing of angular momentum of halos through major mergers. The 
authors define a parameter / which should have on average a 
value of / 0.42 to make their model reproduce the distribu- 
tions found in simulation. However we find a mean value / ~ 1 . 
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which is larger than that required by their model, suggesting an 
inconsistency. 

All presented results are equally valid for major mergers 
(M1/M2 < 4, Ml > M2) and show no dependence on the mass 
ratio or mass of the merging objects supporting a self similar 
build up of structure. The use of the self-consistent initial 
conditions presented here will allow to test consistently the 
merger scenario for the formation of elliptical galaxies in the 
hierarchical universe without covering the whole parameter 
space and provides appropriate statistical weights for different 
initial conditions. 
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